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Summary 

Tests were conducted to evaluate the feasibility 
of using a scale-model experiment situated in an ane- 
choic facility to investigate long-range sound propa- 
gation over ground terrain. At a nominal scale factor 
of 100:1, attenuations along a linear array of sbc mi- 
crophones colinear with a continuous-wave type of 
source were measured over a wavelength range from 
10 to 160 for a nominal test frequency of 10 kHz. 
Grazing-incidence angles ranged from 0° for both 
sound source and microphones set flush to about 15° 
for maximum source and microphone elevations and 
for minimum propagation range. Most. tests were 
made for a hard model surface (plywood), but lim- 
ited tests were also made for a soft model surface 
(plywood with felt). 

For grazing-incidence propagation over the hard 
surface, measured and predicted attenuation trends 
were consistent for microphone locations out to be- 
tween 40 and 80 wavelengths. Beyond 80 wave- 
lengths, however, significant variability was observed 
between test runs. A contributing factor to the run- 
to-run variability at the more remote microphones 
was the relatively long term disturbances in the prop- 
agation medium caused by natural ventilation of the 
anechoic facility. There was also some evidence of 
extraneous propagation-path contributions to data 
irregularities at the more remote microphones. 

An analytical model-sensitivity study at the high- 
est test frequency of 12.5 kHz indicated that same- 
direction systematic errors on the order of ±1 mm 
in both source and microphone elevations could 
generate up to 2.5-dB shifts in “relative” excess at- 
tenuation trends. For equal elevation errors in oppo- 
site directions, however, a cancellation effect appar- 
ently occurred that reduced the sensitivity by about 
50 percent. For the soft surface, no comparable 
sensitivity was found. 

The results of this study suggest that for an ex- 
perimental arrangement consisting of a continuous- 
wave pointlike source, the source and microphone el- 
evations above a hard-surface model should be accu- 
rate to within at least ±0.5 mm for a nominal test 
frequency of 10 kHz. Also, careful attention should 
be given to source design to ensure that extrane- 
ous radiation is at least 40 dB below that from the 
source radiation orifice. If these conditions are sat- 
isfied, then the experiments of the type described in 
this report should prove fruitful for validating ana- 
lytical models for predicting long-range propagation 
over ground terrain, provided that surface impedance 
can be appropriately modeled. 


Introduction 

Reliable estimates of sound attenuation for near- 
grazing-incidence propagation to a surface are 
required for a number of applications. Since high- 
frequency acoustic disturbances suffer the greatest 
atmospheric attenuation, surface-induced losses are 
of importance primarily for lower frequencies. Be- 
cause of the numerous physical effects that must 
be taken into account, theoretical prediction is a 
challenging problem. Over the past two decades, 
numerous investigators have contributed to the de- 
velopment of mathematical models for near-grazing- 
incidence propagation. A thorough review of this 
work is given by Pao et al. in reference 1. The 
analyses by Rudnick and Ingard (refs. 2 and 3, re- 
spectively) provided the first fundamental consider- 
ation of sound radiating from a pointlike source in 
the vicinity of an infinite plane of finite impedance. 
In 1975 Chien and Soroka (ref. 4) extended the anal- 
ysis to include all boundary impedances likely to be 
encountered. Thus, the analytical development for a 
point source radiating above a plane boundary with 
locally reacting impedance appears to be complete at 
the present time. 

Further development to account for the effect of 
surface irregularities on sound propagation over ter- 
rain at low grazing angles has been pursued by sev- 
eral investigators. A summary of this work, given 
by Howe in reference 5, calculated the effect of a 
random distribution of small-scale surface irregulari- 
ties consisting of cylindrical or hemispherical bosses. 
Howe showed that the net effect is the production 
of surface waves that are effective in penetrating 
the ground shadow zone predicted for a smooth sur- 
face of the same “effective” impedance. The effect 
of large-scale surface features with finite impedance 
boundaries such as ridges and hills is currently being 
studied by Pierce et al. (See refs. 6, 7, and 8.) The 
refraction of sound by time-invariant mean velocity 
and sound speed gradients has been considered in 
previous studies (see ref. 9) and in a recent textbook 
(see ref. 10). 

All these mathematical models have been devel- 
oped for idealized terrain geometries, temperature 
and wind gradients, and ground impedances. It is 
practically impossible to perform a parametric val- 
idation of these models on the basis of measure- 
ments conducted in complicated terrain geometries 
and meteorological conditions actually encountered 
outdoors. One way to circumvent this difficulty is 
to develop scale-model experiments that will allow 
controlled, systematic variations of the most impor- 
tant geometric and physical parameters. Such exper- 
iments could validate current simplified models and 



provide the basis for further refinements appropriate 
for full-scale outdoor application. 

In reference 11, Hutchins et al. conducted prop- 
agation experiments on a 1:80 scale model with a 
sound source operating over the frequency range from 
8 to 100 kHz. Hard asphalt surfaces were modeled 
with sheet aluminum. Institutional grass was mod- 
eled by expanded polystyrene that was sanded and 
covered with a single layer of tissue paper. The main 
thrust of this work was to use model results directly 
to predict full-scale excess attenuations. Generally, 
good correlations were observed between model and 
full-scale interference frequencies. 

Other scale-model investigations to determine 
wind and terrain effects on long-range sound prop- 
agation have been conducted in special wind-tunnel 
facilities for studying atmospheric boundary-layer ef- 
fects. One such work, reported by Tachibana and 
Yoshihisa in reference 12, used discrete frequencies 
ranging from 12.5 to 100 kHz to study propagation 
over scale-model distances of about 1.8 m. By heat- 
ing or cooling the model surface they were able to 
separate the effects of thermal and wind gradients 
near the model surface. Their results were qualita- 
tively consistent with those observed in field inves- 
tigations. Another investigation reported by Ander- 
son et al. in reference 13 made use of a spark source 
to study propagation over model surfaces. The 
unique advantage of their facility was control of wind 
velocity and boundary-layer characteristics. 

Two model ground surfaces were examined in 
this experiment. One model, designed to simulate a 
hard asphalt surface, consisted of a 1.6-mm-thick alu- 
minum plate installed on a turntable arrangement in 
the test-section floor of the tunnel. The other model 
surface, which simulated grassland, consisted of 
2.5-cm-thick extruded polystyrene with light tis- 
sue paper glued to the surface. The tunnel test- 
section dimensions were 2.4 m wide by 1.82 m high 
by 18.2 m long. The hard reflecting walls re- 
quired a pulse-type source with appropriate data 
editing to remove unwanted reflections. The pulse 
source was driven by an electrical discharge and 
was characterized by relatively high peak pres- 
sure, run-to-run variability, and a continuous spec- 
trum dominated by frequencies ranging from 10 to 
50 kHz. Although a pulse-type source will im- 
prove signal-to-noise ratios, special care is required 
to account for source nonlinearity, run-to-run vari- 
ability, atmospheric absorption, and transducer di- 
rectionality associated with the high-frequency con- 
tent of the pulse spectrum. The authors of ref- 
erence 13 claim to have adequately accounted for 
these factors in their data acquisition and analysis 
procedure. 


The results of reference 13 were stated as inser- 
tion loss of the model surface (excess attenuation) 
in 1 /3-octave bands and for a fixed microphone lo- 
cation of 2.2 m from the source with angles of graz- 
ing incidence ranging from 0.7° to about 8°. For the 
most part, the theory predicted the general trends for 
propagation over hard and soft flat surfaces with no 
wind. The authors note that the theoretical doubling 
of pressure at a hard plane surface was not consis- 
tently achieved at low frequencies. In the deep min- 
ima associated with destructive interference, mea- 
sured levels were higher than predicted levels by as 
much as 20 dB for the hard surface but by only 4 dB 
for the soft surface. The large discrepancy for the 
hard surface was attributed to microphone directiv- 
ity variations and system noise. Also, theory over- 
predicted the constructive interference maxima by 
as much as 4 dB for both the hard and soft surfaces. 
When a model hill with and without wind was tested, 
the results were far less satisfactory. In this case the 
theory did not predict any details of the measured 
surface insertion losses. 

In summary, the authors of reference 13 suggested 
that the results were encouraging. They further sug- 
gested that some of the discrepancy between theory 
and experiment could be due to bending waves in- 
duced in the aluminum plate and to directionality 
characteristics of the microphones. Clearly, the dis- 
crepancies between theory and experiment for the 
simple flat surface under a homogeneous station- 
ary medium should be thoroughly understood and 
eliminated before proceeding to models that include 
surface irregularities. 

This paper reports the results of an experiment 
designed to gain a better understanding of environ- 
mental control requirements for scale-model simula- 
tion of long-range sound propagation over ground 
terrain at small grazing-incidence angles. Specifi- 
cally, tests were conducted on a 50-m^ expanse of 
model surface located inside an anechoic chamber to 
minimize reflections. A discrete- frequency (10 kHz) 
continuous-wave sound was radiated from a pointlike 
source whose elevation above the model surface could 
range from 0 to about 2 wavelengths. Two different 
model-surface impedance conditions were used in this 
experiment. The baseline surface consisted of high- 
quality marine plywood to simulate a relatively hard 
surface condition. A soft-surface condition was simu- 
lated with a 3-mm-thick felt covering that was spread 
over the plywood surface. An existing data acqui- 
sition system was used to acquire attenuation data 
along a linear array of microphones colinear with 
the source. Synchronous signal averaging was used 
to improve signal-to-noise ratios for the relatively 
small, and therefore insensitive, microphones. The 
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microphone elevations ranged from 0 to about 
1.5 wavelengths above the surface. From these mea- 
surements, attenuation characteristics of the surfaces 
were determined. Comparisons of measured and 
calculated “relative” excess attenuation trends were 
obtained based on estimated values of the model- 
surface reflection factors. Finally, recommendations 
are given based on the experience gained during this 
study to improve scale-model experiments further for 
long-range-propagation model validation. 

Symbols and Abbreviations 

Values are usually given in SI units and, where 
considered useful, also in U.S. Customary Units. 


Ae 

excess attenuation 

(■^e,calc)rj 

calculated excess attenua- 
tion at ri relative to that 
at reference location 

(■^e,meas),.j 

measured excess attenua- 
tion at ri relative to that 
at reference location ri^ref 


free-field source strength 

Ag 

source strength in pres- 
ence of model surface 

A/D 

analog-to-digital 

converter 

c 

sound speed 

D 

range from source to 
microphone receiver 

DVM 

digital voltmeter 

dB/dd 

decibels per distance 
doubling 

erfc 

complex error function 

F{w) 

boundary-loss factor 

FFT 

Fast Fourier Transform 

/ 

frequency 

Hr 

elevation of microphone 
receiver 

Hs 

source elevation 

^rel 

relative humidity, percent 

i 

= v^ 

ki 

free-space (atmosphere) 
wavenumber, u/c 


h 

wavenumber in model 
medium 

Mic 

microphone 

^ref 

reference sound pressure 
level 

(P.)„ 

acoustic pressure at 
field point in presence 
of surface 

■Ptest 

test sound pressure level 

Pb 

barometric pressure, mm 
Hg 

Q 

image source strength 

R 

normal-incidence reflec- 
tion factor 

Rp 

plane wave reflection 
coefficient 

Ri 

model-medium flow 
resistivity 

n 

path length from source 
to microphone receiver 

T2 

path length from image 
to microphone receiver 

SPL 

sound pressure level, dB 

ASPL 

change in sound pressure 
level, dB 

Ta 

ambient temperature, °C 

Ti 

inside temperature, °C 

TTL 

terminated transmission 
line 

V 

outdoor windspeed, knots 

w 

numerical distance 
parameter 

a 

attenuation constant 

0 

normal-incidence 

admittance 

9 

acoustic resistance 


normal-incidence 

impedance 


characteristic impedance 
of model medium 

p 

density of air 

T 

propagation constant of 
model medium 

<t> 

grazing-incidence angle 


X 

'll) 


U) 

Subscripts: 

calc 

ref 


acoustic reactance 

angular locations of 
circular-arc microphones 
relative to radial-array 
microphones, deg 

angular frequency 


calculated 

reference 


Experimental Setup and Procedure 

General Description 

The design of a scale-model, long-range propaga- 
tion experiment requires compromises in the avail- 
able space, ambient environment control, source 
power and radiation characteristics, microphone sen- 
sitivity and size, and data acquisition and analysis 
system. This experiment was exploratory in nature 
and did not permit the design and development of 
the most desirable source and microphone charac- 
teristics. In particular, source power, stability, and 
mechanical isolation from the model surface were 
marginal. Also, microphone sensitivity was com- 
promised for small size to minimize directionality 
effects. It was anticipated at the beginning of the 
design phase that these compromises could be com- 
pensated for to some extent by means of an ex- 
isting sophisticated data acquisition and analysis 
procedure. 

The design and layout of the experimental setup 
is shown in the plan view diagram of figure 1. The 
microphone locations are represented by the open cir- 
cles with distances from the source (shaded symbol) 
shown in wavelengths (A) at the nominal test fre- 
quency of 10 kHz. The experimental objective was to 
simulate long-range sound propagation over a model 
surface between a point source and microphone ar- 
rays at low, audible frequencies. This simulation was 
done by means of a scale-model arrangement that 
trades off a reduction in geometric dimensions for 
an increase in frequency. Thus, k\D remains con- 
stant, where ki is the free-space wavenumber and 
D is the propagation range. To avoid contaminat- 
ing reflections, the model surface was required to be 
effectively infinite and bounded by a nonrefiecting 
enclosure. This was accomplished by installing the 
model surface on the floor of the Langley Anechoic 
Noise Facility and by providing curved shrouds along 
the surface boundaries to minimize edge diffraction. 


Realistic modeling of ground impedance requires 
that the model-surface impedance be equal to that 
of the full-scale surface over the frequency range of 
interest. It is unrealistic to expect available materi- 
als to exhibit the same impedance over the test fre- 
quency range of interest as that of specific ground 
surfaces ranging from vegetation-covered soil to con- 
crete. In this initial experiment, impedance modeling 
was a secondary consideration. The model-surface 
choices were 19-mm-thick (3/4-in.), high-quality ma- 
rine plywood and the same surface covered with 
3-mm-thick (1/8-in.) felt. 

A radial and circular-arc array of microphones 
centered on the source location was positioned as in- 
dicated in figure 1. The radial microphones were 
spaced at successive distance-doubling points with 
the first microphone located about 10 wavelengths 
from the source, based on a nominal test frequency 
of 10 kHz. This frequency choice approximated 
a model/full-scale dimension ratio of 1:100. Five 
circular-arc microphones 120 wavelengths from the 
source were installed nominally 18° apart (symmet- 
rically located with respect to the linear array). The 
purpose of the circular-arc array was to monitor po- 
tentially troublesome behavior attributable to source 
directivity, model-surface vibration, and residual- 
edge diffraction or reflections. To encourage desyn- 
chronization in residual-edge-diffracted wave arrival 
times, the radial array was skewed with respect to 
the surface diagonal (i.e., 23.4° with respect to the 
longer surface edge). Note that microphone 5 serves 
as both a radial and a circular-arc receiver. 

The model surface was assembled from 15 mod- 
ules, each with dimensions 1.2 m by 2.4 m (4 ft by 
8 ft), as indicated in the sketch of figure 2(a). Each 
module consisted of a continuous sheet of marine ply- 
wood attached to a wood support structure around 
its periphery. The modules were supported in the 
anechoic chamber by evenly spaced, steel rails run- 
ning parallel to the concrete floor that normally sup- 
port sections of a “false” plywood floor about 0.3 m 
above the concrete floor of the facility. The total el- 
evation of the model surface was about 0.6 m above 
the concrete floor or 22.2 cm (8.75 in.) above the few 
remaining sections of the false plywood floor. This 
installation allowed sufficient access space for source 
driver, microphone preamplifiers, and cabling. As- 
sembly of the complete model surface was conducted 
by one-at-a-time placement of the modules on the 
support rails. Gaps between adjacent modules were 
filled with epoxy-based filler that was sanded flush 
with the module surfaces. Upon each placement, 
shims were used to achieve an elevation as nearly 
equal as possible to the previously installed mod- 
ules by use of a carpenter’s bubble level. In this 
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manner, point-to-point slope changes, over the length 
of the level (1.22 m), were maintained to within 
±0.3°. Upon completion, the structure could support 
the weight of test personnel, an important consider- 
ation because it was necessary to change the source 
and microphone elevations manually. Lateral com- 
pression forces on the entire assembly were supplied 
by a series of tensioning cables strung through the 
module supports. 

As previously mentioned, the model-surface pe- 
riphery was terminated by a curved shroud to min- 
imize contamination by edge-diffracted waves. The 
details of this structure are illustrated in figure 2(b). 
The critical feature of the shroud design was the 
smooth transition from the fiat surface to the curved 
surface with minimum surface discontinuity. The ra- 
dius of curvature, 20.3 cm (8.0 in.), was large rela- 
tive to the longest wavelength of interest (i.e., 3.8 cm 
at 9 kHz). Past experience with such shrouds, re- 
ported in reference 14, indicated that contamination 
from edge-diffracted waves was minimized by this 
technique. 

The source consisted of a high-frequency audio 
tweeter (5 to 15 kHz) available at audio supply 
outlets. The tweeter was modified and attached to 
a brass tube (with an outside diameter of 1.9 cm 
and an inside diameter of 0.64 cm) by means of 
an adapter section as illustrated in the sketch of 
figure 3. The inside diameter of the brass tube 
was such that only plane wave propagation could be 
supported up to. about 15 kHz, the upper frequency 
limit of the driver. Plane wave propagation through 
the tube ensured cylindrically symmetric radiation 
from the source. The 10-cm length of the brass tube 
allowed a maximum source elevation of 5.8 cm. For 
zero source and microphone elevations, this source 
arrangement was capable of generating 76 to 95 dB 
over the frequency range from 9 to 13 kHz at the 
reference microphone location. Thus, an SPL range 
from about 52 to 71 dB would be expected at the 
most distant microphone location for a hard surface. 

The microphones used were 3-mm-diameter (1/8- 
in.) condenser microphones. The standard adapters 
furnished by the manufacturer to connect the mi- 
crophones to their respective preamplifiers allowed a 
maximum elevation of 3.6 cm above the model sur- 
face without the larger diameter preamplifier section 
protruding into the sound field. In the standard con- 
figuration recommended by the manufacturer, the 
lower limit of the dynamic range of the microphone 
and preamplifier system is 76 dB. The signal level at 
the most distant microphone location was estimated 
to range from 50 to 55 dB for zero source and mi- 
crophone elevations. Time-domain synchronous av- 
eraging was expected to extend the lower limit of the 


dynamic range to provide the additional signal-to- 
noise ratio needed at the more distant microphones. 

To permit the easy removal of both the sound 
source and the microphones, circular sections or 
plugs of plywood were machined to fit snugly into 
the appropriate cutouts of the model surface. (See 
fig. 4.) The undersides of the mounting plugs were 
equipped with a cam and locking arrangement to al- 
low a mechanically hard connection to the model sur- 
face. To remove the assembly for calibration and 
maintenance, a special spanner-type wrench was de- 
signed to disengage the cam mechanism, thereby 
permitting the assemblies to be lifted from their 
seated positions. During test operations, the source 
driver was embedded in a mass of modeling clay 
(about 2 kg (1 lb)) to reduce source driver radia- 
tion into the access space underneath the model sur- 
face. Despite the added mass, measurements with 
the source removed from the model surface indicated 
a decrease of only about 20 dB in the SPL upon plug- 
ging the source radiation orifice. Although this level 
of extraneous radiation from the source driver was 
of some concern, it was anticipated that the model 
surface would shield the “test space” from the un- 
desired radiation since the driver would be located 
underneath the model surface. 

Calibration Apparatus 

Figure 5 is a schematic diagram of the apparatus 
for performing frequency response comparisons of the 
microphone measurement systems. The device con- 
sisted of a 5 mm by 1 mm rectangular waveguide run- 
ning through the center of a rectangular aluminum 
block. The size of the waveguide ensured that only 
plane acoustic waves were transmitted at frequencies 
below about 20 kHz. A reference microphone and a 
measurement microphone were mounted flush on the 
waveguide walls opposite each other. The waveguide 
was then connected to a sound source at one end, 
and the exit end was terminated in a length of tub- 
ing to reduce end reflections inside the channel. In 
this manner, comparison frequency responses of all 
microphone systems could be obtained over a time 
period of about 1 hr. These calibrations, repeated 
at intervals during the tests, were used to evaluate 
the overall repeatability of the data acquistion and 
analysis system. An absolute calibration was per- 
formed for all the microphones using a commercially 
available calibrator operating at 1.0 kHz. 

Data Acquisition and Analysis 

The on-line data acquisition and analysis system 
used in this experiment was an adaptation of that 
developed for the source monitoring system of the 
spinning mode synthesizer in the Langley Aircraft 
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Noise Reduction Laboratory as described in ref- 
erence 15. The instrumentation layout is shown 
schematically in figure 6. For a given test frequency, 
source input current and sound pressure measure- 
ment systems were under direct computer control, 
as indicated by the instruments marked by an as- 
terisk. Upon establishing an acceptable source level, 
the sound field at the microphone locations was sam- 
pled in a systematic manner by switching the multi- 
plexer under computer control. Microphone signals 
were conditioned by the autogain amplifier and 50- 
Hz-bandwidth tracking filter to provide optimum in- 
put to the A/D converter. The maximum sampling 
rate of the A/D converter was 40 kHz which per- 
mitted about 2.7 samples per cycle at a frequency 
of 15 kHz. Source driver stability was monitored by 
performing a measurement at the reference micro- 
phone each time that a measurement was completed 
at a more distant microphone. (That is, sound pres- 
sure levels were measured in the sequence: Mic 1, 
Mic 2, Mic 1, Mic 3, ..., etc.) An FFT analyzer was 
used to monitor broadband spectra from microphone 
system outputs as indicated in the instrumentation 
schematic diagram. 

Time-domain, synchronous averaging reduces the 
contribution of incoherent noise in proportion to 
where N is the number of averages. Thus, 
the low-end dynamic range of a measurement system 
can theoretically be improved indefinitely. However, 
there are practical limits on the time available for 
the process. In this experiment, a sufficient amount 
of synchronous averaging was done to provide an 
effective signal enhancement of about 20 dB. 

To accomplish signal enhancements of up to 
20 dB, 200 blocks of pressure time history data, each 
consisting of 160 data points, were phase synchro- 
nized to a common time origin and then averaged. 
For a nominal test frequency of 10 kHz, this proce- 
dure allowed up to 7200 cycles of the fundamental 
waveform to be processed for each set of 200 data 
blocks. From the stored “averaged digitized wave 
form,” the amplitude and phase of a sinusoidal least 
squares fit were calculated. The stored calibration 
factors were then applied to calculate the SPL and 
phase. The above procedure was repeated five times, 
and then the resulting five SPL’s and phases were av- 
eraged to obtain a final value for a given microphone 
location. 

Because of direct or natural ventilation of the ane- 
choic facility, temperature and relative humidity were 
not controlled. Thus, temperature and relative hu- 
midity tended to correlate with outdoor meteorolog- 
ical conditions. These parameters were recorded by 
sensors located on one wall of the facility. Absolute 
accuracies for the temperature and relative humidity 


readings were not greater than 1°C and 1 percent, 
respectively. 

Analysis 

Propagation Model 

As an aid to evaluating the scale-model propa- 
gation experiment described in this paper, relative 
trends of the measured and computed excess attenu- 
ation are compared. The computations are based on 
the analysis of a point source radiating in the pres- 
ence of a finite impedance surface of infinite spatial 
extent. The simplified formulation given by Ches- 
sell (see ref. 16) was deemed adequate since, for the 
parameter range of interest in this investigation, the 
surface wave is insignificant. 

Figure 7 depicts the relevant source and micro- 
phone geometry for a point source near an infinite 
plane. The upper region represents the atmosphere; 
the lower region, the model medium. The source and 
microphone elevations are taken as Hg and Hr, re- 
spectively; the grazing-incidence angle is <f) and the 
range is D. The distances from the source and image 
to the microphone are denoted by r\ and r 2 , respec- 
tively. The wavenumber of the disturbance is A:i in 
the atmosphere and k 2 in the model medium. At- 
mospheric absorption is neglected but the possibil- 
ity of extended reaction of the model medium is in- 
cluded. For discrete-frequency emission from a point 
source, the pressure {Ps)^ at the microphone con- 
sists of the superposition of the contributions from 
the direct and reflected fields. Thus, the total pres- 
sure at ri, with the temporal dependence suppressed 
for convenience, is given by 

^ exp(ifciri) + ^ exp(iA:ir 2 ) (1) 

Here, As is the source strength and Q is an effective 
image source strength, where 

Q = Rp + F{w) (1 - Rp) (2) 

Note that Q is defined in terms of the plane wave 
reflection coefficient Rp and the so-called boundary- 
loss factor F{w). For surfaces characterized by ex- 
tended reaction, i.e., internal wave motion parallel to 
the surface, the plane wave reflection coefficient for 
a medium with an effectively infinite depth is 

sin 4> — (3\l — {ki/k 2 )‘^ cos^ <^J ^ 

~ r iT/y 

sin <t) + — {k\/k 2 )‘^ cos^ <^J 
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which, for a locally reacting surface, reduces to 


_ sin </) — /? 
sin (3 


( 4 ) 


More elementary analyses replace the image 
source strength Q with the plane wave reflection co- 
efficient Up. However, for grazing incidence, (/> = 0° 
and ri = r 2 - Consequently, Rp = —I regardless of 
the value of /3. Thus, the elementary analyses predict 
the complete absence of an acoustic field for 0 = 0°. 
The boundary-loss factor F{w) is a function of the 
numerical distance w which for an extended reacting 
surface is given by 


2ikiV20^ ki cos^0\ 

{l-Rpf cos2</) V ^2 / 


and for a locally reacting surface becomes 


w = ^ikir 2 


(sin <p + l3)^ 
1 + /3 sin 4> 


( 5 ) 


(6) 


The boundary-loss factor describes phenomena in the 
vicinity of the impedance boundary associated with 
spherical wave propagation at near-grazing-incidence 
angles. When k\r 2 is small or j3 is small (or both), 
F{w) « 1 and Q 1, independent of Rp, thus 
resolving the anomaly that results when Rp is naively 
used in place of Q in equation (1). For the other 
limiting cases, i.e., with high frequencies, longer 
ranges, and larger surface admittances, then it; >> 1 
and F{w) ^ 0. Equation (2) then reduces to the 
familiar plane wave result. Chessell calculates F{w) 
from an asymptotic series good for absolute values of 
w less than 1 or greater than 10 and neatly avoids the 
range of lit;| between 1 and 10 over which the current 
experiment takes place. Therefore, it was necessary 
to evaluate F{w) directly using a result from Chien 
and Soroka (ref. 4) as follows: 


F{w) = 1-1- i\fim exp(— It;) erfc (7) 


For later application it will be convenient to write 
equation (1) with the direct-wave contribution iso- 
lated as follows: 


(P,) = (i+riQ exp|«,(,,-roi| 

^ n [ V2 J 

( 8 ) 

The bulk of current experience suggests that 
acoustic propagation over real soil surfaces can be 
adequately described by analytical models incorpo- 
rating the locally reacting assumption. However, 


for the model surface consisting of felt on plywood, 
significant extended reaction effects may exist. 

A comparison of equations (3) and (4) suggests 
that extended reaction can alter the incidence an- 
gle at which Rp changes sign for certain critical 
combinations of the incidence angle and surface 
admittance. For maximum source and microphone 
elevations, the grazing-incidence angles in this ex- 
periment ranged from about 0.6° at the most dis- 
tant microphone to about 10° for the nearest micro- 
phone. Because of the possible sign change in re- 
flection coefficient at the more distant microphone 
locations, an effort was made to estimate /? as well 
as k 2 at the test frequencies of interest. These es- 
timates, over the test frequency range from 9.0 to 
13.0 kHz, were obtained with the aid of impedance 
measurements over a limited frequency range from 
0.5 to 3.0 kHz. In addition, flow resistance mea- 
surements were obtained for the felt covering. These 
measurements were used with a semiempirical model 
for porous materials developed by Delany and Bazley 
(see ref. 17) to provide the desired estimates of j3 and 
k 2 over the test frequency range. This procedure will 
be discussed further in the “Results and Discussion” 
section. 

Procedure for Correlating Experiment With 

Theory 

Excess attenuation at a point ri in the presence 
of a surface is defined as an insertion loss of the 
surface expressed in decibels relative to the free-field 
pressure at the same point. In the literature, it 
is generally assumed that the source strength and 
directivity are invariant for a change in the radiation 
environment. If the source strength should change 
because of the presence of the surface, the “apparent” 
excess attenuation follows from equation (2) as 



+ 20 log 


1 + 


^Qexp [iki(v2 -ri)]| 

( 9 ) 


where Ay is the free-field source strength. For 
As = Af (i.e., no change in source strength due to the 
surface), the first term is zero and the result would 
reduce to the “true” excess attenuation because of 
ground effects alone. 

In the present experiment, not only did As ^ Af 
but also As changed slowly with time because of 
acoustic driver degradation. It was therefore incon- 
venient to compare experiment and theory directly 
on the basis of equation (9), which implies a knowl- 
edge of the ratio As/Af. Instead, measured and cal- 
culated relative excess attenuations were compared. 
This procedure will now be described. 
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A knowledge of Ag alone would allow the sound 
pressure level at a point ri to be calculated from 
equation (1) as follows: 


(SPLcalc)ri = 20 log 


Pref 


= 20 log 



|1 


+ (’•iA 2)<3 exp[ifei(r 2 - rj)] I 


( 10 ) 


where Pj.g{ is the reference sound pressure level 
(20 /iPa). The source strength Ag can be estimated 
via equation (10) by inserting measured values of 
(SPL)^^ at a reference location j-ef- Thus, a rel- 
ative excess attenuation at any other field point can 
be calculated; i.e., 


(^e,calc),.j - (SPLcale)ri 201 og ^ 

^’"^-Qexp [tfci(r2,ref “ ’’l.ref)] 

( 11 ) 


■20 log 


1-f 


’’2,ref 


Measured relative excess attenuations can be com- 
pared with the calculated values from equation (11) 
by performing the following operation on the mea- 
sured sound pressure levels: 


(>4c,meas)^j — (SPLineas)^ “ 20 log ^ ^ 

'( 12 ) 

It should be noted that both the calculated and mea- 
sured relative excess attenuations as given by equa- 
tions (11) and (12) will be zero at the reference 
location. In the experiment, this condition is actually 
achieved by normalizing the experimental and cal- 
culated values to zero at microphone 1, a procedure 
which effectively accomplishes the same result as for- 
mally solving equation (10) for the source strength 
Ag. Clearly, this procedure allows only measured and 
calculated trends of excess attenuation to be com- 
pared. However, for this particular experiment, the 
difference between calculated and measured relative 
excess attenuations is indicative of inadequacies in 
the analytical model and/or sensitivity to systematic 
errors in measured parameters. 


Results and Discussion 


Comparisons of Microphone Frequency 
Response 

Comparisons of the system frequency response 
which includes source control, data acquisition, and 
microphone measurement systems were made during 
the tests. In figure 8 the results of five representative 


comparison responses, obtained with the specially 
constructed calibration device described previously, 
are presented. The tests were conducted during a 
time period of several days to provide an indication 
of system stability. 

The response comparisons are presented in fig- 
ure 8 for the six microphone systems used in the 
radial-array measurements. The data were obtained 
using the same computer-controlled source and data 
acquisition system used in the scale-model experi- 
mental setup. Responses of each microphone chan- 
nel are presented for five different dates as indicated 
in the figure key. The key also includes ambient 
temperature, relative humidity, and barometric pres- 
sure. The responses were measured at frequency 
increments of 1.0 kHz and ranged from 9.0 to 
13.0 kHz. 

With two exceptions, the maximum variability 
among the six microphone systems, as defined by 
the dashed parallel lines, ranged from 0.6 dB for 
microphones 4 and 5 to 0.8 dB for microphone 1. 
At 12 kHz, the responses for microphones 1 and 2 
ranged up to 1 dB. When such excursions occurred, 
they were scrutinized as possibly being indicative of 
equipment malfunction. 

The conclusion drawn from this exercise is that 
the system response, including microphone system 
and data acquisition system, was invariant with re- 
spect to both time and frequency to within about 
0.8 dB. Thus, 0.8 dB sets a lower limit below which 
the variability associated with sound propagation 
over the model surface may be difficult to resolve. 
Also, 20-Hz-bandwidth spectra of the background 
noise obtained with an FFT analyzer, as seen by 
the microphones installed on the model surface, indi- 
cated levels from 52 to 54 dB over the test frequency 
range. For a perfectly reflecting surface, a pure tone 
level of 80 dB would be required at microphone 1 to 
produce a level of 56 dB at microphone 6, the most 
distant microphone. Synchronous signal averaging 
was expected to provide at least 20-dB improve- 
ment in the signal-to-noise ratio at the most distant 
microphone. 

Measured Model-Surface Impedance 

The normal-incidence impedances or, equiva- 
lently, the complex reflection factors (see eq. (5) 
with = 0°), for the plywood and felt-covered ply- 
wood test specimens were measured. Photographs 
of two square test specimens (^^=^5 cm by ^5 cm) are 
shown in figure 9. These measurements were per- 
formed in an impedance tube with the same cross- 
sectional dimensions as the test specimens. The spec- 
imen impedance was determined from the standing 
wave pattern generated by normally incident and 
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reflected waves at the test specimen surface. For 
highly reflecting surfaces with reflection factors above 
about 0.95 (standing wave ratios of 32 dB or greater), 
direct measurement of the standing wave parameters 
is prone to imprecision and inaccuracy because of 
ambient noise interference with null level measure- 
ments. A multipoint sampling technique has been 
developed at the Langley Research Center that cir- 
cumvents the need for direct measurements in sharp 
nulls as required by the standing wave method. The 
technique essentially allows a reconstruction of the 
standing wave pattern from a selected number of 
pressure and phase measurements. The impedance 
is calculated from the reflection factor obtained from 
the best fit of a one-dimensional propagation model 
to the data. Comparisons between reflection factor 
magnitudes using this method and the standing wave 
method for highly reflective materials such as steel 
indicate that the multipoint method can consistently 
measure reflection factor magnitudes in excess of 0.95 
with a repeatability of ±0.005. 

Figure 10 shows the measured normal-incidence 
impedances (normalized by pc) plotted against fre- 
quency for both model surfaces. (See the figure key 
for symbol definitions.) For comparison, the mea- 
sured impedance trend for a soil surface as reported 
by Embleton et al. in reference 18 is represented by 
the short-long dashed curves. As might be expected, 
both impedance components for the hard and soft 
model surfaces decrease in magnitude with increasing 
frequency. The soil-surface impedance components 
have a similar frequency dependence. However, the 
impedance components for the plywood surface de- 
crease at precipitous rates compared with those for 
the soil surface, at least up to the maximum measure- 
ment frequency of 3.0 kHz. The resistive component, 
in particular, appears to decrease below the soil resis- 
tance mean value at 3.0 kHz. When the felt covering 
is added, the impedance components are in better 
agreement with those of the soil surface, especially 
from about 2.0 to 3.0 kHz. 

The challenge of devising a model-surface imped- 
ance to appropriately scale a “real” surface at a cho- 
sen scale factor is dramatically illustrated by these 
data. For example, the indicated graphical extrap- 
olations of the plywood-surface data (short dashed 
curves) suggest that at a scale factor of about 20, the 
reactance component may be a fair approximation of 
that of soil at 0.5 kHz as indicated by the shaded 
areas centered on 0.5 and 10 kHz; however, the re- 
sistance component is not well-modeled. In a later 
section, these impedance data will be used in con- 
junction with a semiempirical model to extrapolate 
the model-surface impedance into the test frequency 
range of interest. 


Acoustic Propagation Test Data 

In this section, acoustic propagation test results 
will be presented and discussed. During the course of 
the experiment, 296 runs were obtained over a time 
period of about 2^ months. The test configurations 
consisted of source and microphones all flush, all at 
maximum elevation, all at one-half elevation, and 
various combinations of source or microphones flush 
and the other at one-half or maximum elevation. Test 
frequencies ranged from 9.0 to 13.0 kHz in 0.5-kHz 
increments. Most of the runs were conducted for 
the hard surface with several repeat runs throughout 
the time period to document test data variability. 
All the basic data are presented in two formats. 
Along the radial microphone array, the distribution 
of sound pressure levels in decibels is plotted versus 
log(fciD), where D is the range from the source to 
a particular microphone. Along the radial array, the 
microphone at location 1 (microphone 1) served as a 
reference. At the nominal test frequency of 10 kHz, 
D ranged out to about 160 wavelengths. Along the 
circular-arc array, the distribution of sound pressure 
in decibels relative to microphone 5, jointly shared by 
the radial and circular-arc arrays, is plotted versus 
angular distance from the radial array. 

Run-to-run, random variations of the measured 
SPL’s were observed among the circular-arc micro- 
phones and the more distant radial microphones. 
These variations suggested a combination of extrane- 
ous propagation paths and temporal changes in the 
propagation medium. A more thorough discussion of 
these factors will be reserved until after presentation 
of representative data. 

To illustrate the data quality, the results of a 
repetitive sequence of six runs are shown in figures 11 
and 12 for frequencies from 10.0 to 12.5 kKz con- 
ducted on 2 different days over a 4-hr timespan and 
for both source and microphones flush with the hard 
surface. The “a” and “b” parts of each figure show 
the radial and circular-arc SPL distributions, respec- 
tively. On the radial distribution plots, 6-dB/dd 
reference lines have been superimposed to indicate 
expected spatial attenuation rate due to spherical 
spreading for an infinite, perfectly reflecting surface 
and nonabsorbing medium. In figure 11(a), which 
shows the radial SPL distributions, there are several 
features worthy of note. First, source output changes 
by about 20 dB as the frequency is incremented in 
0.5-kHz steps from 10.0 to 12.5 kHz. This is caused 
mostly by the change in loading on the driver because 
of resonant behavior of the source adapter. Second, 
at test frequencies of 12.0 and 12.5 kHz for which the 
reference levels at microphone 1 exceeded 90 dB, the 
attenuation rate is measurably greater than 6 dB/dd, 
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in a consistent pattern, at the more distant micro- 
phones. Third, the other test frequencies, whose ini- 
tial reference levels were less than about 85 dB, tend 
to exhibit attenuation rates greater than 6 dB/dd, 
but with significant random variation both within 
a given run and among runs for the more distant 
microphones. 

Figure 11(b) shows the SPL variation along the 
circular-arc array. The pressure levels in this fig- 
ure are plotted in decibels relative to that at mi- 
crophone 5, which is at the center of the circular- 
arc array and is shared with the linear array. (See 
fig. 1.) Angular variation in the sound pressure lev- 
els for a given test frequency is seen to be great- 
est between the —36° and —18° microphone lo- 
cations, reaching nearly 12 dB for run 1 at the 
10-kHz test frequency. With the exception of run 1, 
the variations for the two positive angular locations 
were 3 to 5 dB less than those for negative angular 
locations. 

Figure 12 shows a repeat of the test sequence 
shown in figure 11 but conducted 4 days later. In 
figure 12(a) the most obvious change compared with 
figure 11(a) is the 5-dB decrease in reference level of 
the 12.0-kHz test frequency and the 2-dB increase 
for the 11.0-kHz test frequency. For the remain- 
ing test frequencies, the reference levels are approx- 
imately the same as those for the morning runs of 
figure 11(a). There are, however, changes in lev- 
els from 3 to 5 dB at the more remote microphone 
locations. 

A comparison of figure 12(b) with figure 11(b) re- 
veals that the circular-arc variation in sound pressure 
levels changed significantly between the test dates. 
At the —36° angular location, all levels changed rel- 
ative to microphone 5 by up to 14 dB. At other angu- 
lar locations, especially the positive ones, the changes 
were not nearly as great. In figure 12(b) the variabil- 
ity at the —36° location was greater than that at 
the same angular location on the earlier test date, 
whereas at the other circular-arc locations the vari- 
ability was generally less than that on the first test 
date. Figure 12(b) also suggests that the run-to-run 
variability increased with angular distance from the 
radial array. 

In the remainder of this report, attention will be 
concentrated on two representative frequencies, 9.0 
and 12.5 kHz, near the extremes of the frequency 
range investigated. The first set of three figures 
will show “raw” data for the 9.0-kHz test frequency 
and for several repeat runs. Two source and micro- 
phone configurations will be discussed. The first 
configuration had the source and microphones all 
at zero elevation above the model surface, and 
the second configuration had the source at its 


maximum elevation of 5.84 cm and the microphones 
all at their maximum elevation of 3.56 cm. Most 
of the data will be for the hard-surface condition. 
Limited soft-surface results will be presented for 
the elevated source and microphones. When the 
source and microphones were at zero elevation on 
the soft surface, large attenuation rates precluded 
adequate signal-to-noise ratios for meaningful re- 
sults. Following the 9.0-kHz data, the same se- 
quence for the 12.5-kHz data will be presented. Al- 
though the circular-arc data were collected for all 
runs, no further presentations of these data will 
be made because they exhibited the same variabil- 
ity discussed previously throughout the experiment. 
Finally, the results will be cast in relative excess 
attenuation form and correlated with calculations 
based on the model-surface impedance extrapola- 
tions discussed previously. 

Figure 13 shows hard-surface radial SPL distri- 
butions for four repeat runs at a test frequency of 
9.0 kHz and for the source and microphones at zero 
elevation. The 10-dB change in levels between runs 
226 and 229 is due to the replacement of the source 
driver between those runs. Note in the figure key 
that these data were taken over a time period of 
10 days during which the outdoor windspeed varied 
by 4 knots, the indoor temperature varied by 2°C, 
and the indoor relative humidity varied by 9 percent. 
As was the case for the previous radial distribution 
data, the levels generally attenuate at a rate greater 
than 6 dB/dd in a consistent manner through micro- 
phone 4. Again, the deviations at microphones 5 and 
6 are typical of the data variability observed at the 
more distant microphones shown in previous figures. 

Figure 14 shows the SPL distribution data similar 
in format to that of figure 13 except for six repeat 
runs with the source and microphones at maximum 
elevations. For this configuration, an interference 
maximum occurs between microphones 1 and 2. Note 
that these data were taken on the same day and 
that the temperature varied by about 1°C, typical 
of morning-to-afternoon changes. The radial SPL’s 
are seen to vary by not more than 1 dB at the 
various microphone locations, with the exception 
of microphone 6 which exhibits a 3-dB variation. 
In general, there was less run-to-run variation for 
elevated source and microphones. 

The behavior of radial SPL distributions for 
the source and microphones at maximum elevations 
above the soft surface is illustrated in figure 15 for 
seven repeat runs. Here, four out of seven runs show 
relatively good radial SPL distributions. These four 
runs (270, 273, 276, and 279) were all conducted on 
the same day ranging from the morning to afternoon 
during which an air temperature change of about 1°C 
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was recorded. It is of interest to note that run 285, 
the last run of the series and conducted on the same 
day as the other four runs listed above, is in good 
agreement except at microphones 4, 5, and 6. Also, 
the outdoor windspeed had been recorded at 13 knots 
during the previous hour even though it had subsided 
to 8 knots during the actual run. 

Figures 16 to 18 show radial SPL distributions 
for the 12.5-kHz test frequency in the same sequence 
and format as in the previous figures for the 9.0-kHz 
test frequency. These data confirm the same general 
trends observed for the 9.0-kHz case. Again there 
was large variability in the circular-arc SPL distribu- 
tions, which are not shown. As with the 9.0-kHz case, 
the radial distributions for the hard surface tend to 
exhibit increased variability at the more remote mi- 
crophone locations with the source and microphone 
elevations at zero. For the elevated source and micro- 
phone configuration, run-to-run data consistency at 
the remote microphone locations improves. The ad- 
dition of a soft surface does not significantly change 
the data variability for the maximum source and 
microphone elevation. 

Data Variability and Possible Causes 

Possible causes of the data variability observed in 
figures 1 1 to 18 will now be discussed. The variability 
can be classified into two categories. First, there was 
a relatively slow systematic decrease in the output 
level of the source with time for constant driver 
voltage (which was always set at the maximum “safe” 
level). Second, there was random variability observed 
for the circular-arc array and to a lesser extent at 
microphones 5 and 6 along the radial array. This 
random variability observed at microphones 5 and 6 
was manifested by deviations from the 6-dB/dd trend 
at the more distant microphones and the variation at 
these same microphones in repeat tests. 

The time scale for changes in driver output, at- 
tributable to driver degradation and to temperature- 
dependent resonance behavior, was such that the 
variable driver output was not a problem over a 
time period of 20 min required for a typical run. 
Driver degradation over operating time periods of 
8 to 16 hr usually mandated driver replacement. 
The criterion for driver replacement was failure to 
maintain a level of at least 76 dB at microphone 1 
for a test frequency of interest and for zero source 
and microphone elevations. This criterion was based 
on an assumed, minimum level of 52 dB at micro- 
phone 6. The ambient broadband noise level ranged 
from 52 to 54 dB as measured by an FFT ana- 
lyzer set for a 20-Hz bandwidth and connected to 
the microphone multiplexer output. (See fig. 6.) 
Synchronous time-domain signal averaging extended 


the effective system noise floor down to about 20 dB 
on a unit bandwidth basis. 

The random variability was likely associated with 
time-varying environmental influences on the direct, 
reflected, as well as extraneous, propagation paths. 
The time scale for these variations ranged from ap- 
proximately 1 hr to several days. Thus contamina- 
tion from competing propagation paths other than 
those emanating from the source radiation orifice is 
believed chiefly responsible for observed data vari- 
ability. This possibility appears to be supported by 
the results of some diagnostic tests conducted with 
power applied to the source but with the radiating 
orifice plugged to reveal “leakage” radiation. Com- 
parisons of level changes at microphones 1 and 6 
for unplugged and plugged conditions at 12.5 and 
13.0 kHz indicated only a 7-dB level change at mi- 
crophone 6 for a 24-dB level change at microphone 1. 
The level changes at microphone 6 were not limited 
to the broadband noise floor because a spectral peak, 
with the source orifice plugged, was clearly evident 
above the background noise floor even without sig- 
nal enhancement. Inspections conducted at the edge 
shroud termination near the false facility floor re- 
vealed detectable sound propagating in the access 
space between the model surface and anechoic cham- 
ber floor. All these tests were made after a quick 
remedy was implemented to reduce leakage radiation 
by embedding the source driver in about 2 kg (1 lb) 
of modeling clay. 

Assuming that such extraneous propagation 
paths alluded to above were present, then it would 
appear possible that thermal gradients induced by 
convection currents or by the normal thermal cy- 
cling of the building structure could alter the relative 
propagation speeds and thus cause time-dependent 
interference with the desired signals. In addition, 
time-dependent refraction was possibly induced by 
thermal gradients near the model surface. Time- 
dependent refraction can cause signal fluctuations 
at the distant microphones even in the absence of 
extraneous propagation paths. Also, residual-edge- 
diffracted and reflected waves from small hard objects 
in the anechoic chamber such as light bulbs and elec- 
trical conduits may have been contributors to spu- 
rious propagation paths. However, a simple scatter- 
ing calculation for a worst-case situation, consisting 
of a 25-mm-diameter cylinder placed 1 m from mi- 
crophone 6, indicated a scattered amplitude at least 
20 dB below the incident wave. 

To recapitulate, the temporal data variability at 
the most distant microphones was apparently due 
to the combined effects of source “leakage” into 
undesired propagation paths coupled with time- 
varying convective and refractive effects, driven by 
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thermal cycling and natural ventilation of the build- 
ing structure. Marginal signal strength relative to 
the noise floor of the broadband instrumentation sys- 
tem wets not perceived to have been a contributing 
factor. Although these effects significantly degraded 
the data quality at the more distant microphone lo- 
cations, the data from the less distant microphones 
(i.e., locations 1 through 4) were relatively stable 
and consistent. Therefore, the overall results were 
encouraging. 

Comparison of Theory With Experiment 

In this section, comparisons between experimen- 
tal and analytical results, as discussed previously in 
the “Analysis” section, will be presented. Represen- 
tative results will be discussed for 9.0 and 12.5 kHz, 
for zero and maximum elevations of source and mi- 
crophones, and for the hard and soft model surfaces. 
To help relate the relative excess attenuation results 
to the radial, absolute SPL’s already discussed, inter- 
mediate results will be presented in which all levels 
are referenced to microphone 1. Also, these plots 
include calculated attenuations based on both hard- 
and soft-surface reflection factors. Finally, the same 
results will be shown in terms of measured and calcu- 
lated relative excess attenuations for various choices 
of the surface flow resistivity parameter. As indi- 
cated in the “Analysis” section, the boundary-loss 
factor F(ty) is a key parameter that describes the 
spherical wave interaction with the finite impedance 
boundary condition. In the test frequency range 
of this experiment, the values of surface impedance 
needed to evaluate F(tt;) were extrapolated by a pro- 
cedure to be described below. 

Because it was not possible to measure the model- 
surface impedance with the available equipment in 
the test frequency range of interest, it was necessary 
to extrapolate the measured impedance to these fre- 
quencies. To accomplish this extrapolation in a ra- 
tional manner, an attempt was made to fit the mea- 
sured data by a semiempirical model produced by 
Delany and Barley. (See ref. 17.) This model per- 
mits the acoustic properties of a porous medium to 
be estimated based on measured values of material 
flow resistivity. For completeness, the equations rep- 
resenting the model are given as follows: 


= l + 0.0571(p//Ei)“°’^^'‘ 

(13) 

X = -0.087(p//fli)-0-^32 

(14) 

: = ki [ 0 . 189 ( p /// 2 i )-®-® 9 ^] 

(15) 

: A:i [l + 0.0978(p//i?i)~° '^‘“] 

(16) 


Equations (13) and (14) give the normalized charac- 
teristic resistance 0 and reactance x in terms of the 
flow resistivity parameter Equations (15) and 

(16) give the attenuation constant a and wavenumber 
k 2 in terms of the wavenumber for air k\ and the flow 
resistivity parameter. The characteristic impedance 
and propagation constant of the material then be- 
come, respectively. 


= ^ + 

(17) 

T = a + ik2 

(18) 


The normal-incidence impedance of a bulk-reacting 
material (felt in the present case) of depth I on a 
highly reflecting surface (plywood) was adapted from 
a well-known formula for a transmission line termi- 
nated in an impedance characterized by a reflection 
factor i?. (See eq. (1.04) in ref. 19 by Zwikker and 
Kosten.) Making use of the estimated values of 
and r from equations (17) and (18), respectively, the 
surface impedance is 

[(1 + fi) cosh(r/) + ^c(l - 7?) sinh(r/)] 
^[(l-hi2) sinh(r/) + ^c(l ““ ^) cosh(r/)] 

(19) 

from which the surface admittance /? is obtained by 

^=l/e (20) 

Values of (3 calculated from equations (19) and (20) 
and values of k 2 taken from equation (16) were used 
in equation (4) to estimate i?p. Then, Rp and k 2 were 
used in equation (5) for estimating w. Then, F{w) 
was calculated from equation (7) and, finally, Q from 
equation (2). For an extended reacting surface, equa- 
tions (3) and (5) are used instead of equations (4) 
and (6) to calculate Rp and respectively. Directly 
measured values of |i?p|, obtained in an impedance 
tube at low frequencies, were compared with calcu- 
lated values at higher frequencies for both hard and 
soft model surfaces. 

Equations (13) to (20) were applied to the ply- 
wood model surface treated as a single layer of in- 
finite depth and to the felt-covered plywood treated 
as a two-layer medium. The two-layer medium was 
modeled with the aid of the modified, terminated 
transmission line (TTL) model as specified by equa- 
tion (16). A measurement of flow resistance for the 
felt was obtained by standard means, the results of 
which are presented in figure 19 for two felt test spec- 
imens. The median value for two samples at zero 
velocity was taken to be as indicated in figure 19. 
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The plywood terminating surface in the two-layer 
medium was characterized as being either a perfectly 
reflecting surface or a reflection factor taken from 
the extrapolation described above for the uncovered 
plywood surface. 

Figure 20 shows the extrapolated reflection factor 
magnitude into the test frequency range of interest 
for the hard and soft model surfaces. For the scope of 
this study, only the reflection factor magnitude was 
extrapolated. The normal-incidence reflection factor, 
as opposed to the impedance, is preferred for this ex- 
trapolation because the separate effects of reflection 
factor magnitude and phase changes have greater 
physical meaning than a change in the real and imag- 
inary impedance components. Reflection factor mag- 
nitude relates to the energy reflected by the surface 
and therefore to the interference null depths, whereas 
the phase part tends to control the interference null 
locations. Because it is unrealistic to expect a scale- 
model surface to match both impedance components 
to those of the full-scale surface, it is perhaps more 
meaningful to make a choice between reflection factor 
magnitude and phase than between resistance and 
reactance. 

In figure 20 the short-long dashed curve shows 
the trend of the measured reflection factor magni- 
tude for a soil surface as calculated from the data 
of reference 18. The real and imaginary parts of 
the measured impedance are presented in figure 10. 
This curve is included in figure 20 for reference pur- 
poses only. No attempt will be made to match the 
soil-surface and model-surface reflection factor mag- 
nitudes in this experiment. Extrapolations of re- 
flection factor magnitudes for the hard surface were 
obtained from the semiempirical model (ref. 17) by 
Delany and Bazley (eqs. (13) and (14)) applied to 
the hard surface treated as a porous medium of in- 
finite depth. The plot shows the results for three 
different values of the flow resistivity parameter. Ap- 
parently, the measured reflection factors, represented 
by the square symbols, are best described by a flow 
resistivity between 20 x 10® and 40 x 10® kg/m®-sec. 

Also in figure 20, calculated reflection factor mag- 
nitudes are presented for the soft-surface specimen 
for two different flow resistivities of the felt. The 
felt-covering thickness was taken as 3.1 mm (1/8 in.). 
The reflection factor magnitude for the soft surface 
did not differ significantly when the plywood reflec- 
tion factor was changed from 1.0 to the extrapolated 
values from figure 20. No reasonable value of the flow 
resistivity, including the possibility of a plywood- 
surface reflectivity less than 1, allowed the TTL 
model to describe the measured data (diamond sym- 
bols). The TTL model curves are relatively close 
below 4 kHz but diverge above 4 kHz for the two 


indicated resistivities. The higher flow resistivity is 
associated with the lower reflection factor in the high- 
frequency range. 

Another surprising feature of the measured reflec- 
tion factor magnitude for the hard- and soft-surface 
specimens is the occurrence of more scatter for the 
soft-surface reflection factor magnitude. Generally, 
more scatter would be expected for the higher reflec- 
tivity associated with the hard surface. These results 
indicate that the multipoint sampling method is ca- 
pable of producing repeatable results at reflectivities 
approaching 1.0. Thus, it appears that the scatter or 
variability in the data for the soft-surface specimen is 
somehow associated with the felt-on-surface compos- 
ite arrangement. One potential source of variability 
for such test specimen configurations consisting of 
contiguous layers of different materials in intimate, 
but unbonded, mechanical contact is the likelihood 
of small air spaces or nonuniform mechanical con- 
tact at the layer-to-layer interface. (See ref. 20 by 
Smith and Parrott.) It should be noted that the 
soft-surface model does not provide an acceptable 
“scaled impedance” for any known outdoor surface 
and therefore serves only as a lower limiting value 
for this scale-model validation experiment. 

Figure 21 shows the results of a parametric eval- 
uation of the boundary-loss factor F{w) for a range 
of the numerical distance parameter w (see eqs. (6) 
and (7)) likely to be encountered. Only log|F(tw)| 
versus |t/;| is shown here for brevity. Referring to 
equation (2) for the image source strength Q, it will 
be recalled that when F{w) approaches 1 (i.e., when 
the surface impedance is high or kiR is small), the 
image strength factor is approximately -|-1 and is in- 
dependent of the value of Rp. For large k^R and 
smaller surface impedances, F{w) approaches zero 
and the plane wave result is retrieved. In figure 21 
the range of numerical distance for the hard and soft 
model surfaces is indicated for both flush (zero el- 
evation) and maximum source and microphone el- 
evations. Clearly, this experiment did not explore 
the entire range of numerical distance; however, rep- 
resentative segments of the low, middle, and high 
ranges of the numerical distance magnitude were in- 
vestigated. Also note that increased source and mi- 
crophone elevation increases the |u;| range for the 
hard surface, whereas for the soft surface the ele- 
vation has no significant effect on the range of |u;|. 
In figure 21(b) relatively large oscillations occur in 
F{w) for values of |t/;| ranging between 5 and 20. 

Figure 22 shows the radial-array attenuation data 
of figure 13 for 9.0 kHz at zero source and micro- 
phone elevations. The attenuation is normalized to 
microphone 1 which, in effect, removes run-to-run 
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variability from the data. Also on the plot, calculated 
attenuation curves have been superimposed for a 
perfectly reflecting surface (i.e., (5 = 1; see eq. (2)) 
and for high and low values of flow resistivity ob- 
tained from the extrapolations for the plywood sur- 
face as shown in figure 20, For this particular case, 
calculated results for the high value of flow resistivity 
and a perfectly reflecting surface are essentially the 
same. For the low value of flow resistivity, the calcu- 
lated attenuation gradually deviates from the curve 
for a perfectly reflecting surface to a value about 
1.5 dB greater at the most distant microphone. With 
the exception of microphone 5, the high flow resis- 
tivity value correlates with the upper bound of the 
measured attenuations. The roll-off of the calculated 
attenuation for the low value of flow resistivity tends 
to correlate with the tendency of the measured at- 
tenuations to roll off from the 6-dB/dd curve at the 
more remote microphones. 

Figure 23 shows the data of figure 22 expressed 
in the format of relative excess attenuations. The 
sign convention for excess attenuation chosen here is 
the same as that used by Chessell in reference 16. 
Therefore, positive excess attenuation indicates that 
a sound pressure level at a given microphone is lower 
than it would be in the presence of a nonabsorbing 
surface. Essentially, this data presentation format 
removes the 6 dB/dd and thereby allows the “non- 
spherical spreading effects” to be studied separately. 
Thus, in the presence of a perfectly reflecting rigid 
surface bounded by an inviscid medium, the rela- 
tive excess attenuation for this configuration would 
be zero which, from the comparison with the almost 
perfectly reflecting surface in the previous figure, is 
approximated by the solid line labeled with a flow 
resistivity of 80 x 10^ kg/m^-sec. Again, the cal- 
culated relative excess attenuation for a surface flow 
resistivity of 20 x 10^ kg/m^-sec has also been in- 
cluded. Note that the term “relative excess attenu- 
ation” means that the calculated and measured val- 
ues have been forced to agree at microphone 1. This 
effect produces the same results that would be ob- 
tained if the measured SPL at microphone 1 were 
used to calculate the source strength. Differences 
between calculated and measured values at other mi- 
crophone locations indicate inadequacies of the theo- 
retical model, sensitivity to systematic measurement 
errors, or random scatter due to extraneous noise. 

Figure 23 suggests that no reasonable value of 
surface flow resistivity allows the predictive model 
to describe the behavior of the measured relative 
excess attenuation in an average sense. Despite the 
measured data variability, the general trend of the 
data suggests significant absorption at zero grazing- 
incidence angle. 


Figures 24 and 25 present the same plot sequence 
as that in figures 22 and 23 but with the source 
and microphone elevations at their maximum values. 
Figure 24 again shows no change in the predicted 
attenuation between a perfectly reflecting surface 
and the highest value of plywood flow resistivity of 
80 X 10® kg/m^ -sec. Note that the data variability is 
much less when the source and microphones are ele- 
vated. Also, even though the data trend is accurately 
predicted, the theory is 2 to 3 dB below the data 
trend. Again, it should be emphasized that perfect 
agreement at the reference location is forced to occur, 
which effectively provides the source strength As- 
This procedure causes any positioning error sensitiv- 
ity at microphone 1 to manifest itself as an “appar- 
ent” systematic discrepancy between the calculated 
and measured results at other locations. This consti- 
tutes a disadvantage of this method of comparison. 
Nevertheless, in the absence of all analytical model 
inadequacies and systematic measurement errors, the 
analytical and measured results would agree. 

Figure 25 shows the same data in the relative ex- 
cess attenuation format and with calculated curves 
for two additional surface flow resistivities as indi- 
cated in the figure. Apparently, the data trend cor- 
relates slightly better with the intermediate value of 
flow resistivity. Note also that the data for the el- 
evated source and microphones are on the negative 
side of the calculated curve in contrast with the zero- 
elevation source and microphones, in which case the 
data tended to be on the positive side of the cal- 
culated curve. These shifts of the calculated curves 
with respect to the data again suggest sensitivity to 
measured parameters used as input to the analytical 
model. As an example of such sensitivity, the effect 
of small changes in source and microphone elevations 
on the calculated relative excess attenuation will be 
discussed next. 

Figure 26 shows the effect on the calculated re- 
sults when source and microphone elevations are var- 
ied from their nominal values of 5.84 and 3.56 cm, 
respectively. For same-direction elevation changes of 
about 2.5 mm, the curves shift by 1.5 to 2.5 dB, In 
fact, positive elevation increments of 2,5 mm for both 
the source and the microphone bring the predicted 
behavior into excellent agreement with run 241. El- 
evation errors for the source and microphones for 
the hard surface were not greater than about 2 mm. 
However, small, randomly distributed changes in sur- 
face slope up to about ±0.3° due to the surface in- 
stallation procedure could have added significantly to 
the relative source and microphone elevation error. 

Figure 27 shows a comparison to measured and 
calculated attenuations for the maximum source and 
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microphone elevations above the soft-surface model. 
As can be observed, there is a dramatic increase 
in the attenuation as compared with that for the 
hard surface. (See fig. 24.) Data for zero elevation 
above the soft surface were not meaningful because 
of the rapid attenuation of the acoustic level below 
the instrumentation plus ambient noise floor. It is 
of interest to note that the measured attenuations in 
figure 27 apparently confirm a conservative dynamic 
range of at least 50 dB in that several attenuations in 
the neighborhood of 50 dB relative to microphone 1 
are obtained. Again, these observations seem to 
suggest that the instrumentation noise floor may not 
have been the cause of random variability, but rather 
contributions from spurious propagation paths. The 
bracketing values of 30 x 10^ and 65 x 10^ kg/m^- 
sec for the soft-surface flow resistivity (see figs. 19 
and 20) generate essentially the same prediction for 
the attenuation profile. Interestingly, the general 
trend of the measured attenuations appears to be 
well-predicted. 

Figure 28 shows the same data plotted in the rel- 
ative excess attenuation format. Although not in- 
cluded in this figure, small changes in the source and 
microphone elevations had no discernible effect on 
the calculated results for the soft surface. Again, ex- 
cept for the large data variability, the calculated re- 
sult appears to describe the average behavior of the 
data. These results are interesting because they sug- 
gest that the agreement between the math model and 
experiment tends to be less sensitive to microphone- 
positioning error than is the case for the hard-surface 
condition, even in the presence of severe contaminat- 
ing influences. 

The final set of five figures (figs. 29 to 33) 
shows comparisons of measured and predicted rel- 
ative sound pressure level distributions and relative 
excess attenuations along the radial microphone ar- 
ray for the test frequency of 12.5 kHz. This frequency 
is near the upper limit of the frequency range inves- 
tigated. These figures are presented in the same for- 
mat as that for the 9.0-kHz data. The results are 
similar to those for 9.0 kHz and tend to support 
the same conclusions; consequently, the discussion of 
these data will be abbreviated. 

Radieil sound pressures in decibels relative to 
microphone 1 for six repeat runs at zero source 
and microphone elevations, taken from figure 16, are 
shown in figure 29. Superimposed on the data are 
calculated levels for surfaces characterized by flow 
resistivities of 20 x 10® and 80 x 10® kg/m®-sec 
and by a perfectly reflecting surface (Q = 1) as 
indicated in the figure. As in the 9.0-kHz case, the 
flow resistivity of 80 x 10® kg/m®-sec is effectively a 
perfectly reflecting surface. 


As previously noted for the 9.0-kHz test fre- 
quency, the data repeatability is reasonably good 
except for microphone locations 4, 5, and 6 which 
again exhibit greater variability. Referring to 
figure 16, three of the runs (227, 230, and 233) were 
made in the afternoon of the same day but during 
a time of modest change in outdoor windspeed and 
temperature. Runs 230 and 233 are in good agree- 
ment for all microphone locations except for num- 
ber 4, for which there is a difference of about 2 dB. 
Starting at microphone 3, run 227 deviates most 
widely from the otW runs. The remaining runs (207, 
224, and 254) were made on different days, and for 
run 207 the air temperature was about 3°C to 4°C 
below that for the other runs. As previously noted, 
there is no problem with the measurement system 
noise floor for these data because the lowest abso- 
lute level is 60 dB, or over 30 dB above the effec- 
tive measurement system noise floor. These observa- 
tions provide further support for the possibility that 
spurious propagation paths were responsible for the 
variability. 

In figure 30 the same data are presented in rel- 
ative excess attenuation format. Clearly, data vari- 
ability tends to mask any discernible trends; how- 
ever, the bulk of the data points fall in the positive 
region of relative excess attenuation. Thus, as in the 
corresponding configuration for the 9.0-kHz case (see 
fig. 23), there is a suggestion that the surface is less 
than a perfectly reflecting one at grazing incidence. 

For the source and microphones at maximum ele- 
vations, the attenuations normalized to reference mi- 
crophone 1 for the hard-surface condition are shown 
in figure 31. Note that the reference microphone is 
now well within the first interference null. The out- 
lying data points fall 2 to 5 dB above the nearest 
calculated curve. Also, note that at this frequency a 
surface flow resistivity of 80 x 10® kg/m®-sec gener- 
ates a result that is quite different from the perfectly 
reflecting surface {Q = 1) and seems to be in better 
agreement with the data than either the perfectly re- 
flecting surface or the flow resistivity surface of 20 x 
10® kg/m® -sec. This statement is confirmed by fig- 
ure 32 which shows these same data in the relative 
excess attenuation format. As with the correspond- 
ing plot for the 9.0-kHz data (fig. 25), the data fall 
on the negative side of the calculated curves. 

Figure 33 illustrates the effect of same-direction 
increments of 0.64 mm (0.023 wavelengths) in source 
and microphone elevations from the nominal, mea- 
sured values of 5.842 cm and 3.556 cm for the source 
and microphones, respectively. This elevation-error 
magnitude is well within the range expected for this 
experiment. The figure clearly demonstrates that 
changes in relative excess attenuation of 2 to 3 dB 
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can result for same-direction systematic errors on 
the order of 1 mm at this particular test frequency. 
However, if the errors are in the opposite direction, 
i.e., source error positive and microphone error neg- 
ative, calculations (not shown in the figure) indicate 
that some cancellation occurs, thus reducing the ef- 
fect by about 50 percent. The results of figures 26 
and 33 suggest that for hard-surface scale models, op- 
erating at nominal test frequencies of 12.5 kHz, dif- 
ferences between measured and predicted relative ex- 
cess attenuation of 1 dB can easily result from source 
and microphone elevation errors of 0.5 mm. 

Because of possible thermal gradients near the 
model surface that may arise from thermal cycling 
of the building structure, a calculation was made to 
determine the relative effects of same-direction ele- 
vation errors and a 0.5°C thermal gradient extend- 
ing 2 cm above the model surface. For maximum 
source and microphone elevations, such a thermal 
gradient was found to cause negligible phase shift 
compared with that due to the same-direction eleva- 
tion errors of 0.64 mm at a propagation distance of 
40 wavelengths. 

Concluding Remarks 

Tests were conducted to evaluate a scale-model 
experiment for studying long-range sound propaga- 
tion at small grazing-incidence angles ranging from 
0® to 15^^. Both hard- and soft-surface models were 
tested. The hard surface was achieved with high- 
quality marine plywood, and the soft surface was 
achieved by covering the plywood with a layer of 
3-mm-thick felt. For a nominal scale-model operating 
frequency of 10 kHz, a discrete-frequency pointlike 
source radiated sound over an array of microphones 
arranged colinearly with the source. Source and 
microphone elevations ranged within about 2 wave- 
lengths of the model surface. The entire appara- 
tus was located in an anechoic facility. The micro- 
phones were positioned at distance-doubling incre- 
ments from the source, starting at 10 wavelengths 
and going out to 160 wavelengths. 

Numerous repeat runs were performed over a 
time period of about 2 months to document data 
consistency for a range of frequencies and combina- 
tions of source and microphone elevations. Gener- 
ally, the data consistency from run to run was good 
over a propagation range between 40 and 80 wave- 
lengths. Beyond 80 wavelengths, significant variabil- 
ity was evident. Also, variability between runs was 
generally greater when the source and microphones 
were at zero elevation, i.e., flush with the model sur- 
face. Contributors to the variability may have been 
residual-edge diffraction, that is, reflections from ob- 
jects such as light bulbs and time-varying refraction 


in the propagating medium associated with natural 
ventilation of the anechoic facility. In addition, extra- 
neous propagation paths arising from the imperfectly 
isolated source driver located underneath the model 
surface may have aggravated the variability problem. 

General trends of the measured “relative” excess 
attenuations were in good agreement with predicted 
trends for the less-distant microphones. Even when 
significant data variability was present for the more 
distant microphones, the average trends were pre- 
dicted by the analytical model. The most critical 
experimental parameters, as suggested by an ana- 
lytical model-sensitivity study, were same-direction 
source and microphone elevation errors. Such errors, 
which could easily range up to 2.5 mm in this experi- 
ment, were shown to result in up to 3-dB shifts in the 
predicted relative excess attenuation curves for the 
hard-surface model. When source and microphone 
elevation errors were in opposite directions, however, 
a cancellation occurred that reduced the shift by at 
least 50 percent. For the soft-surface model, such 
shifts were not predicted. 

The results of this study suggest that for sound 
radiating from a pointlike source over hard-surface 
models, source and microphone elevation accuracies 
should be better than ±0.5 mm to limit systematic 
discrepancies between measured and predicted rela- 
tive attenuation to about 1 dB for a nominal test 
frequency of 12.5 kHz. Also, careful attention should 
be given to source design to ensure that extrane- 
ous radiation is at least 40 dB below that from the 
source radiation orifice. If these conditions are sat- 
isfied, then the experiments of the type described in 
this report should prove fruitful for validating ana- 
lytical models for predicting long-range propagation 
over ground terrain, provided that surface impedance 
can be appropriately modeled. 


NASA Langley Research Center 
Hampton, Virginia 23665-5225 
September 2, 1987 
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Figure 4. Photograph of front view of sound source and microphone assemblies. 




Figure 5. Schematic diagram of apparatus for calibrating microphone systems (relative frequency responses). 
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Figure 6. Schematic diagram of instrumentation layout. 
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Figure 7. Sketch showing geometry of math model. 
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Figure 8. Frequency response comparisons for microphone systems. 
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(a) Soft surface. (b) Hard surface. 

Figure 9. Photograph of model-surface test specimens for impedance measurements. 
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Figure 10. Comparison of normal-incidence impedances of model surfaces with that of soil surface. 
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(b) Circular-arc SPL distributions. 
Figure 12. Concluded. 
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microphones at elevations of 5.84 and 3.56 cm. 



L distributions along radial array at frequency of 12.5 kHz 
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Figure 19. Measured flow resistance of felt covering. R\ Ri 5.9 X 10^ kg/: 




Figure 20. Comparison of measured and calculated reflection factor magnitudes for hard and soft model 
surfaces. 
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Figure 21. Log of boundary-loss factor versus magnitude of numerical distance. 



Figure 21. Concluded. 








<D 

> 

o 

XJ 



o 

oo 


00 

CM 


CO 

CM 


CM 


CM 

CM 


O 

CM 


00 



L 


O 

i—H 

I 


Lf) 

r-H 

I 


CO 


CM 

I 


u 


00 

-o 


(O 


s: ^ 

< o 

4-> 


CD 

S- 


Figure 24. Comparison of measured and calculated changes in radial SPL distributions at frequency of 9.0 kHz 
for hard surface with source and microphones at elevations of 5.84 and 3.56 cm. 
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Figure 25. Comparison of measured and calculated relative excess attenuation for three values of R\ at 
frequency of 9.0 kHz and inside temperature of 25°C. Hard surface; source and microphones at elevations 
of 5.84 and 3.56 cm. 
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Figure 27. Comparison of measured and calculated changes in radial SPL distributions at frequency of 9.0 kHz 
and inside temperature of 26°C. Soft surface; source and microphones at elevations of 5.84 and 3.56 cm. 





Figure 29. Comparison of measured and calculated changes in radial SPL distributions at frequency of 
12.5 kHz and inside temperature of 25°. Hard surface; source and microphones at zero elevation. 
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Figure 30. Comparison of measured and calculated relative excess attenuation for two values of R\ at frequency 
of 12.5 kHz for hard surface with source and microphones at zero elevation. 
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Figure 31. Comparison of measured and calculated changes in radial SPL distributions at frequency of 
12.5 kHz for hard surface with source and microphones at elevations of 5.84 and 3.56 cm. 


239 


CM UD CX) «-H 
^ ^ LO 
CM CM CM CM 


DO<^0 


04j 


□ O <5d 


a 4> 


CM C7> 
O 


Figure 32. Comparison of measured and calculated relative excess attenuation for two values of i?i at frequency 
of 12.5 kHz for hard surface with source and microphones at elevations of 5.84 and 3.56 cm. 
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Figure 33. Comparison of measured and calculated relative excess attenuation for three source and microphone 
elevations at frequency of 12.5 kHz for hard surface. 
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